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Abstract 

! ! The electroweak baryogenesis is studied in the charge transport 

' mechanism with the vector-like quark model. Introducing an extra 

vector-like up-type quark and a singlet Higgs scalar with the mass of 
, the order of a few hundred GeV, the baryon number generation from 

\Q [ the bubble wall is estimated. We show that this scenario is consistent 

with the measurement of the present baryon to entropy ratio of our 
. universe, if the parameters are in the right region. 

^ : 

Jh ; 1 Introduction 

> : 

^ ■ The baryogenesis scenario in the electroweak scale was suggested in order 

to solve a number of problems in the GUT scale baryogenesis. Within the 
minimal standard model, however, this electroweak baryogenesis scenario 
faces other serious problems. First, CP-violation source coming only from 
the CKM phase is too small to explain the observed baryon to entropy ratio. 
Second, the electroweak phase transition is of the weak first order one in the 
minimal standard model, although the strong first order phase transition is 
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required for the baryogenesis to be realized. Third, after the phase transition 
finishes, the difficult wash out problem exists. In order to solve these serious 
problems, we have to go beyond the standard model. 

Recently, Nelson, Kaplan and Cohen proposed the electroweak baryoge- 
nesis scenario called the charge transport mechanism JIJ. If the electroweak 
phase transition is of the first order, bubbles of the broken phase appear in 
the symmetric phase during the development of the phase transition. The 
charge transport mechanism is the scenario in which quarks reflect from the 
expanding bubble wall in a manner that produces net hypercharge flux in 
the symmetric phase. This hypercharge production is guaranteed by the CP- 
violating interactions inside the wall. Then the produced hypercharge flux is 
converted into the baryon asymmetry through the sphaleron transition active 
in the unbroken phase outside the bubble wall. 

In this paper, we study the electroweak baryogenesis with the vector-like 
quark model. This is a minimal extension of the standard model in which 
extra SU(2) singlet vector-like quarks and one extra singlet Higgs scalar are 
added to the standard model. Historically this model was first introduced 
by Bento, Branco and Parada as one of the simplest spontaneous CP- 
violation models. The vector-like quark model is very attractive to the elec- 
troweak baryogenesis, since it can not only accommodate an extra CP viola- 
tion source, but also makes the first order phase transition of the electroweak 
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theory much stronger. The baryogenesis problem with the vector-like quark 
model was studied by J. McDonald |7j using the spontaneous baryogenesis 
mechanism which is different from the charge transport mechanism studied 
in this paper. 

2 The vector-like quark model 

We consider here an extension of the standard model by adding extra SU(2) 
singlet quarks Ul and Ur as the fourth generation as well as an extra singlet 
Higgs scalar S. These extra singlet quarks are called vector-like, since they 
have left-right symmetric weak isospin and their couplings with the gauge 
fields become vector-like. 

The field content of this model is 

(u d)\ , u% d R , U L , 0, S, 
1 = 1,2,3, a = l, 2,3,4 

where i and a are the generation indices, and <fi denotes the standard Higgs 
doublet. We assume that the additional singlet Higgs S is a complex scalar 
and all the new fields introduced, Ul, Ur = u R and S, are odd under a Z 2 
symmetry, whereas the ordinary fields are even. 

Under these conditions, the SU(2) x U(l) x Z 2 invariant Yukawa couplings 
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are given by 

L Y = -V2(ud) l L {h^d R + f ij( f>u R ) - fiU L U R - V2(f i4 S + f' t4 S*)U L u R + h.c. , (1) 

with additional coupling constants /i, f i4 and f' i4 (i — 1 ~ 3). 

We consider a scenario in which the extra singlet scalar, S, obtains a 
vacuum expectation value at the higher energy scale before the electroweak 
phase transition starts, and later the doublet Higgs, 0, obtains the vacuum 
expectation value during the electroweak phase transition. Therefore after 
the electroweak phase transition, S and are able to have the following 
vacuum expectation values, 

<*> = ^(.°) • {S) = t/°- (2) 

The most general renomalizable SU(2) x U(l) x Z 2 symmetric scalar 
potential is as follows, 

V = V* + V S + V^s, 

V s = S*S( ai + b^S) + (S 2 + S* 2 )(a 2 + b 2 S*S) + b 3 (S* + S* 4 ), 
V4,s = f [ Cl (,S 2 + S* 2 ) + c 2 S*S] , (3) 

where all the coefficients are assumed to be real. With the given profile of 
the expectation values, v and V, expectation value of the phase a can be 
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obtained by minimizing the potential with respect to a, that is, 

(4) 

From this eq.(f|) we understand that the value of a can be nonvanishing 
even for v = before the electroweak phase transition starts. Further- 
more it can be position-dependent inside the bubble wall, during the elec- 
troweak phase transition is undertaking, reflecting the space-dependency of 
v there. The position-dependent complex phase is the manifestation of the 
spontaneous CP violation in this model. There is, however, a problem in- 
herent to the spontaneous CP violation, that is, there are 2 types among 
the 4 independent solutions satisfying eq.(|]), having different CP proper- 
ties. That is, if the phase a + is one of the solutions, the phase a_ satisfying 
sin(a_) = — sin(a + ) is also the solution. Therefore, generation of the bub- 
bles having these different CP phases may cancel the net hyperchage pro- 
duction from the bubble walls. In order to avoid this difficulty, an explicit 
but tiny CP breaking may be helpful. Following 0, if we replace the factor 
(S 2 + S* 2 ) in the potential V^s to (e ia °S 2 + e~ iao S* 2 ), then the difference 
of the potential AV between two solutions of opposite signs, a + and ct_ is 
roughly AV ~ 0(V 4 ) sin(a )| sin(2a±)|. Let us consider (— ) bubbles having 
less potential energy than the surrounding (+) bubbles by AV. In order for 
the (— ) bubbles can grow up to fill the whole universe, the critical radius r c 
of the (— ) bubble in the surrounding (+) bubbles should be less than the 
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cos 



a 2 + b 2 V 2 /2 + Cl v 2 /2 
hV 2 



radius of the universe r# : 

r c <r H . (5) 

The critical radius r c is the ratio of the surface tension a ~ 0(V 3 ) over the 
AV, r c ~ a/AV ~ 1 / (V sm(a )\ sin(2a±)|). The radius of the universe, r H , 
is related to the Hubble constant at temperature T as follows, 

r H = H(T)-\ 

where 

H(T) ~ 20T 2 /M pl . 

Therefore eq.(j^) leads to 

20T 2 

sino;o| sin2a ± | > — — 1(T 16 , (6) 

V M p i 

where we have chosen T = 100 GeV, V = 500 GeV and M pi = 10 19 GeV. 
Therefore the tiny explicit CP breaking of O(10~ 16 ) is enough to keep the (— ) 
bubbles having the same CP properties, eliminating the other (+) bubbles. 

The shape of the wall v(z) can be obtained from the equation of motion 
for 0. Here, we assume that a wall profile can be determined dominantly 
from the potential V$ as 

/ 1 + tanh(z/5„,)\ 
v(z) = v , (7) 



where we consider only the dependency on the coordinate z, normal to the 
bubble wall, but we can get the more general solutions by boosting the wall 
profile ([7p in the x-y plane. Here 5 W represents the thickness of the wall. 

Now, let us consider the Yukawa couplings of the up-type quarks. As the 
weak eigen states, we take the bases in which mass matrix of the down-type 
quarks is diagonalized. In these bases, mass matrix of the up-type quarks 
can be written as 

»<% < s > 

where 

mj = fijV, M 4j = (faS + f' Aj S*). 

This mass matrix can be diagonalized by the bi-unitary transformation, that 
is, by using the two unitary matrices Ul and Ur we have 

UlMU R = M d =(** (9) 

where m = (m u , m c , m t ) and M is the vector-like quark mass Mjj. We assume 
that the vector-like quark and the singlet Higgs are heavier than the standard 
quarks and doublet Higgs. 

3 Computation of the hypercharge flux 

For simplicity, we consider here only top quark and the vector-like quark in 
the computation of the hypercharge flux. What we have to do is to solve 
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the Dirac equation under the position-dependent wall profile and to find 
the reflection coefficients of the quarks from the wall. The two-generations' 
up-quark mass matrix is 

M = ( FVe i«Jl Z F ' Ve -ia(z) (10) 

where / = / 33 , F = / 43 and F' = f' 43 . Using this ^-dependent mass matrix, 
we will write the Dirac equation in the chiral bases. For the stationary state 
(ip(t,z) = ip(z)e lEt ), the Dirac equation reads 



.d_(Mz) \ = _( E Mt\ f^ lR (z) 

dz\^ L {z)) \M -E)\4> 3L (z) 
i 0_(1>*L{z)\__(E MWfc(z)\ 

dz\rhR{z))~ iMt -EJU 2R (z)J [LL) 

where i[>l(r)(z) has two components, corresponding to the top and the vector- 
like quarks. 

We have solved this Dirac equation numerically as a system of ordi- 
nary differential equations with respect to z, without any approximations. 
The reflection coefficient from the right-handed j-quark to the left-handed %- 
quark, R^(E), and from the left-handed j-quark to the right-handed i-quark, 
R lj (E), are given by 

RHE) = (12) 
R ij (E) = ( 13 ) 



where i and j represent the generations, E is the incident energy and z 
denotes the position of the boundary between the wall and the symmetric 
phase. Using these reflection coefficients, the hypercharge flux can be written 

as 

roo roo / | \ 

fy = L dkL L kTdkT {4^) 

E [(Ay)^|i2y(fc L )| 2 + (AY^lR^t 
i 

(rm 3 - f\E%) , (i4) 



hi 

X 



with 



f s {E)i = 
f\E) 



yJtf L + k 2 T + Mf (1 + exp(E s /T)) ' 



yjk\ + k 2 + Mf (1 + exp(E b /T)) ' 
E s = 7 (jk\ + k 2 T + Mf - v w k^ , 

E b = 7 (y/k? L + k 2 T + Mf + v w k L ^j , (15) 

where (AYj^j represents the change of the hypercharge from j-th quark to 
the i-th quark, 7 is the 7-factor 1/ ^1 — v \ corresponding to the wall velocity 
v w , kr and ki are transverse and longitudinal components of the momentum 
in the wall rest frame, and E s (E b ) and M s (M b ) are, respectively, the energy 
and the diagonalized mass of quark in the thermal frame of the symmetric 
(broken) phase. The integral over k T can be done analytically. Then we have 



(Ay)^|i2y(fc L )| 2 + (AY)^ t \R tj (k L ) 



X 



log (l + exp [--y(k L - v wy Jkl - Mf)/T 
— log ( 1 + exp 



l(k L + v w Jkl-Mf)/T 



(16) 



We have to give the phase transition temperature and the masses of the 
quarks in the broken phase for the calculation of the hypercharge flux. We 
set the phase transition temperature T = 100 GeV, and the mass of the top 
quark m t = 174 GeV. The wall width is considered roughly the order of the 
phase transition temperature, <5~ x ~ T ~ 100 GeV, but we have estimated 
here the hypercharge flux for the various choices of the wall width. The 
experimental constraints for the mass of the vector-like quark is not so strict 
[§]. We have tried the cases of M v = 500 GeV and 300 GeV. Most of the 
parameters in the scalar potential V in eq. (3) are outside of the experimental 
verification. Therefore we can choose freely a and b with a restriction of 
a + b = 1 and determine the position-dependent a as cos 2a = a + b(v(z)/v ) 2 
. The restriction above means the disappearance of the spontaneous CP 
violation phase after the electroweak phase transition ends. We also introduce 
rricom and h defined by m com = (F + F')V and h = F' / F . Then, in the case 
of M v = 500 GeV, m com = 300 GeV, and v w = 0.5, we have f Y = O(10~ 6 — 7 ) 
for the various choices of (b, h) = (1.0 ~ 2.0, 0.01 ~ 0.5). The value b = 2.0 
corresponds to the case that the complex phase moves from zero to 7r/2 inside 
the wall. If the value h becomes smaller, the imaginary part becomes larger 
than the real part in the off-diagonal term. The value h = 1.0 corresponds 
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to the real off-diagonal term. 

If we change the off-diagonal value, m com , of the mass matrix, we have 
obtained the following results for M v = 300 GeV, v w = 0.1, h = 0.01, and 
b = 2.0 : fy = O(10~ 8 ~~ 6 ) corresponding to the values of m com = 0.0 ~ 1.2. 
(see Fig.l) 

Fig. 2 shows the wall width dependence of the hypercharge flux in the case 
of Mu = 300GeV and My = 500GeV. We can see that the hypercharge flux 
takes the larger value for the thinner wall. 

4 Generation of the baryon asymmetry from 
the hypercharge flux 

Following the method of Nelson, Kaplan and Cohen M, we compute the 
partial derivative of the free energy with respect to the baryon number. If 
the time scale of the baryon number violating sphaleron transition is much 
longer than that of the weak interactions between reflected quarks, 10~ 26 [s], 
we will be able to think the process of the sphaleron transition is not in the 
thermal equilibrium within the time scale of the weak interactions. Then, the 
baryon number and the lepton number will conserve during the time scale. 
In that situation we can introduce the chemical potentials corresponding to 
hypercharge Y/2, charge Q, baryon number B, B — L (L : lepton number), 
and B = {B± + B^)/2 — (B2 + B\)/2, where Bi is the baryon number of 
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the i-th generation. We have neglected the mixing between the lighter (1st 
and 2nd) generations and the heavier (3rd and 4th) generations, but this is a 
quite reasonable assumption. Using these chemical potentials, we can obtain 
the relation between the chemical potential of the baryon number /ib and 
the hypercharge density py/2 : 

T 2 



-2(24/5 + -)(8(3 + 1) + (24/5 + l)(16/3 + 1) 
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Pb 

pY/2 ~ WTi 

- J 4»bT\ (17) 
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where 



] fOO <j>2 

P=T^ dx (18) 

4tH Jo i + cos h J x 2 + M 2 



One can easily see that the (3 is small enough to be neglected as long as we 
set the mass of the vector-like quark Mjj at a few hundred GeV . Then the 
net baryon number production is 



PB = T ~fjdt 



r.,„h f , . dF 

dB 

9T h r z / v w 

= 13^7-00 dt PY/2(z - V w t) , (19) 

where py/2(x) denotes the hypercharge density at the position which is lo- 
cated at the distance x from the bubble wall. T sph = k(awT) 4 is the 
sphaleron transition rate of the symmetric phase, where aw = g 2 /4:7i, and k 
is a numerical factor which has been estimated by J.Ambj0rn et al. about 
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.1 ~ 1.0 || . The integral over t can be estimated approximately by the 



transport time of the reflected quarks tt as follows: 



/. 



z/v. 




oo 



dtpYji{z - v w i) 



dxp Y/2 (x) 



— oo 



Iytt 



(20) 



v.. 



w 



and the baryon to entropy ratio becomes 



9 ka^T f Y r T 



s 




13 s v. 



(fv/T 3 ). 



(21) 



Then to explain the present baryon to entropy ratio, we need fy ~ 
Iq-5 s i nce typical transport time of the weak scale is roughly 

O(10 1 ~ 3 T _1 ). The solutions of fy ~ 10 -5 6 T 3 have been found in the 
charge transport scenario with the vector-like quark model in the last section, 
so that the scenario can be a good candidate of the model of the baryogenesis. 

5 Conclusion 

In this paper, the charge transport scenario of the baryogenesis has been 
studied using the vector-like quark model. We have adopted the model in 
which extra left-right symmetric, SU(2) singlet up-type quarks Ul and Ur, 
and a SU{2) singlet Higgs scalar S are added to the standard model. The 
possible Higgs potential for the ordinary doublet Higgs and the singlet S 



13 



admits the spontaneous CP violation in which the phase a of S takes the 
non-vanishing vacuum expectation value. During the first order phase tran- 
sition bubbles of the broken phase are generated, and the phase a becomes 
position dependent inside the bubble walls. A tiny but explicit CP violation 
should be, however, incorporated in order to avoid the difficulty of having 
two types bubbles. In the presence of this explicit CP breaking, one type 
of the bubbles creating the baryon number through the charge transport 
mechanism, remains, whereas the other type of the bubbles decreasing it, 
disappears. 

By solving the Dirac equation numerically without any approximations, 
we have estimated the rate of the hyperchage production from the bubble 
walls, in which the presence of the position-dependent CP phase in the mass 
mixing between the top quark and the extra up- type one is important. The 
evaluated hypercharge flux from the wall fy depends on the various unknown 
parameters. However, with the top quark mass of 174 GeV, the vector-like 
quark mass of 500 GeV or 300 GeV, and the phase transition temperature of 
100 GeV, we have solutions giving f Y ~ 10~ 6 for the value of the wall width 
5 W = (lOOGeV)^ 1 . It is the reasonable value for reproducing the baryon 
to entropy ratio of 10~ 10 11 . For the transition of the hypercharge to the 
baryon number, the usual discussion has been extended so as to incorpo- 
rate the extra vector-like up-type quark and the singlet Higgs, in which the 



14 



thermal equilibrium by the weak interactions is violated by the sphaleron 
transition. 

The phase transition dynamics when singlet Higgs is introduced has not 
been given here. Some of the results can be found by using the 1-loop ap- 
proximation and the high temperature expansion Q on how the first order 
phase transition becomes stronger in our case 0, but the more detailed 
study should be required. In the real development of the first order phase 
transition, fusion effect of the bubbles including the temporal change of the 
expanding velocity of the bubbles, should be taken into account ||. This 
kind of analysis is difficult, but very important. Also the detailed analysis 
should be given on the disappearance of the unwanted bubbles having the 
different CP properties. In the analysis of estimating the final value for the 
baryon to entropy ratio, we need to know the exact value of the sphaleron 
transition rate, or the coefficient k in eq.(21), as well as the transport time 
for the various quarks. Improvement on the estimation of these parameters 
are also necessary. 
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Figure Captions 

1. Fig.l: The wall width dependence of the hypercharge flux in the case 
of M v = 500 GeV and m com = 300 GeV for b = 2.0, h = 0.01, v w = 0.1 
and 5 W = T _1 (solid line). The case of My = 300 GeV and m com = 100 
GeV is also shown (dashed line). 

2. Fig. 2: The hypercharge flux as a function of m com = 0.0 ~ 1.2 GeV for 
M v = 300 GeV. Other parameters are the same as in Fig.l. 
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